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Abstract 

We develop a theory of dark matter based on a previously proposed picture, in which a complex 
vacuum scalar field makes the universe a superfluid, with the energy density of the superfluid giving 



y^ \ rise to dark energy, and variations from vacuum density giving rise to dark matter. We formulate a 



nonlinear Klein-Gordon equation to describe the superfluid, treating galaxies as external sources. 
We study the response of the superfiuid to the galaxies, in particular, the emergence of the dark- 



en 

"^! matter galactic halo, contortions during galaxy collisions, and the creation of vortices due to galactic 



rotation. 
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I. INTRODUCTION AND SUMMARY 



The standard model of particle theory postulates the existence of at least one complex 
vacuum scalar field (the Higgs field), in order to generate certain masses in a gauge-invariant 
manner, mainly in the weak sector, through spontaneous symmetry breaking. Direct ex- 
perimental support of this field comes from the discovery of the associated field quantum, 
the Higgs boson [1]. Extension of the standard model to supersymmetric or grand unified 
theories will bring in still more complex vacuum fields. We may therefore suppose that the 
universe is permeated with complex scalar fields, whose manifestation in the large has not 
been the preoccupation of particle theory, but is the main focus of this paper. 
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As we learn from condensed matter physics, a complex scalar field on the macroscopic 
scale gives rise to superfluidity. The field serves as an order parameter that expresses quan- 
tum phase coherence over macroscopic distances, and the superfluid velocity corresponds 
to the gradient of the phase of the complex field. From this point of view, we must view 
the universe as superfluid. It has been proposed [2|, |3|] that this cosmic superfluidity offers 
explanations to both "dark energy" and "dark matter" : 

• Dark energy is the energy density of the cosmic superfluid, which drives an accelerated 
expansion of the universe; 

• Dark matter is the manifestation of local density variations of the superfluid, which 
are detectable through gravitational lensing. 

In this paper, we formulate and apply a phenomenological theory of dark matter based on 
this picture. 
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In Ref.[2l], the dark-energy aspect is explored in a model using a self- interacting scalar 
field coupled to gravity through Einstein's equation with Robertson- Walker metric, in the era 
immediately following the big bang. The Robertson- Walker metric sets the only length scale 
in the problem, and the field potential must change with this scale through quantum field- 
theoretic renormalization. For mathematical consistency, one requires that the potential be 
asymptotically free (i.e., vanish in the small-distance limit,) and this determines it to be 
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the Halpern-Huang potential J4j derived from renormalization-group studies. It is a non- 
polynomial potential with exponential behavior for large fields, a feature leading to the 
generation of an effective cosmological constant that decrease in time like a power, and 
this gives dark energy (accelerated expansion of the universe), while avoiding the usual 
"fine-tuning" problem. 

The same model gives a new scenario of the inflation era [3|]. After the big bang, matter 
(modeled as a perfect fluid coupled to the scalar field) was created in the quantum turbu- 
lence of the cosmic superfluid, in the form of a vortex tangle. It grows and decays in about 
10~^^s, and is described phenomenologically by Vinen's equation originally used in liquid he- 
lium. The inflation era corresponds to the lifetime of quantum turbulence, during which the 
radius of the universe increases by a factor 10^^, and all the matter needed for subsequently 
nucleosynthesis was created. At the end of the inflation era, the usefulness of this model is 



ended, for density variations in the universe became important. The model passes control 
to the standard hot big bang theory, with an important condition: the universe remains a 
superfluid, and all astrophysical processes take place in this superfluid. 

The foregoing describes how the cosmic superfluid comes into being, and its role during 
the inflation era. Now we have to describe it post inflation, in which conditions have become 
somewhat different. 

The big bang era was governed by the Planck scale of lO^^GeV, which was built into 
Einstein's equation through the Robertson- Walker metric. With the emergence of matter 
comes the nucleon mass scale of 1 GeV. This scale ultimately arises in scale-invariant QCD 
through spontaneous breaking of chiral symmetry — an effect referred to as "dimensional 
transmutation". (The observable mass of the universe comes mainly from protons and 
neutrons, and has little to do with the Higgs.) Since the matter length scale is much smaller 
than the Planck scale, it should prevail in the dynamics of stars and galaxies. The decoupling 
between the nuclear scale and the Planck scale is illustrated in Ref. 3|]. In the present work, 
which is largely concerned with the present universe, we simply ignore the influence of the 
Planck scale, and use a simple 0^ theory with constant phenomenological coupling constants. 

In principle, one should use an updated Halpern-Huang potential, with needed renormal- 
ization to bridge the Planck era and the present. Everything that happens between then and 
now — the expansion of the universe, emergence of the matter scale, interactions between 
matter and the scalar field — will contribute to this "running" of the potential, with the 
length scale expanding by 18 orders of magnitude. We bypass this problem here by adopting 
the phenomenological 0^ potential. 

The order parameter originates in Landau's [5|] theory of phase transitions, in which a 
high-temperature phase with less "order" (i.e., higher degree of symmetry) changes into 
a low-temperature phase with more "order" (less symmetry), through what is now called 
spontaneous symmetry breaking. The order parameter is something that gives a measure of 
"order" , and it vanishes above the critical temperature. The generic example is the magnetic 
moment density in a paramagnetic phase transition, which can be represented by a scalar 
field with a symmetric field potential V{(j)), which has only one minimum with = at 
high temperatures, but exhibits two minima at nonzero values = ±0o below the critical 
temperature. This is a phenomenological way to describe spontaneous symmetry breaking. 
In the Ginsburg-Landau theory of superconductivity ^, the order parameter is a complex 



scalar field coupled to the electromagnetic field, and it successfully describes the Meissner 
effect, in which the photon acquires mass inside a superconductor. (This is the first instance 
of the so-called "Higgs mechanism".) We now know that this order parameter corresponds 
to the complex wave function of Cooper pairs in the microscopic BCS theory [7], and the 
Ginsburg-Landau theory has been derived [81] from BCS. However, it continues to serve as 
an indispensable practical tool. 

The Landau approach is also important for understanding superfluidity in liquid helium 
and cold-trapped atoms [9] . The order parameter in this case is the complex wave function 
of the Bose-Einstein condensate, and it obeys the Gross-Pitaevskii equation, an instance of 



the NLSE (nonlinear Schrodinger equation) first used in nonlinear optics 10|. 

In our investigation here, we consider a generic complex scalar field in curved space-time, 
which may correspond to some vacuum field in particle theories, such as supersymmetry or 
grand unified theories. We assume that the field obey the nonlinear Klein-Gordon equation 
(NLKG) in curved space-time. We shall discuss the definition of superfluid density and 
superfiuid velocity in a relativistic setting. Matter will be treated as a perfect fluid with 
given density current, interacting with the scalar field via current-current interaction. In 
this work we have not considered the dynamics of matter itself. 

The NLKG describes a pure superfluid at absolute zero temperature. At higher tem- 
peratures, there will appear a normal fluid, consisting of phonons, field quanta, as well as 



particles coupled to it, existing with the superfluid in thermal equilibrium. The so-ca^ 



led 



WIMPs (weakly interacting massive particles) searched for in particle-physics setups [Uj, 
would be a component of our normal fluid, if they exist, and this would explain how they 
came to permeate all space. However, we have not taken the normal fluid into account in 
this paper. 

The most persuasive evidence for dark matter is the detection of galactic halos via grav- 
itational lensing. In our picture, the dark-matter halo is made up of extra superfluid drawn 
around a galaxy, due to attractions between galaxy and superfluid. The halo can follow the 
motion of a galaxy like a "soliton". It can follow the rotation of a galaxy by developing 
quantized vortices. When two galaxies collide, their halos will merge and undergo contor- 
tions. All these phenomena are part of the superfluid hydrodynamics that follows from the 
NLKG. Observed dark matter formations like the so-called "bullet cluster" [l3J and "train 
wreck" [ij] are, in our picture, manifestations of superfluid hydrodynamics. 



The main topics investigated in this paper are 

• fitting model parameters to cosmological data, 

• dark matter in galactic halos, 

• quantized superfiuid vortices, 

• numerical simulations of dark-matter halos, galaxy collisions, vortices. 

II. NONLINEAR KLEIN-GORDON EQUATION (NLKG) 

We use units in which h = c = 1, and a metric corresponding to (—1,1,1,1) in fiat 
space-time. The vacuum complex scalar field is denoted by 

(x) = F (x) e^"(^) (1) 

The classical Lagrangian density in the absence of galactic matter is 

C = -(7^^9^0*9,0 - V (2) 

where g'^" is the metric tensor, and V (0*0) is the self-interaction potential. The action is 

S = - f dt d^xy/^ {g^''d^(t)*d^(t) + V) (3) 

where g = det((7^iy). This leads to the equation of motion 

^ :5, (v^^-^9.0) -^ = (4) 



^/^ (90* 

We use a phenomenological 0'^ potential, which is the simplest way to maintain a nonzero 
vacuum field Fq: 

V = ^m'-F^Y + Vo (5) 

where Vq is the vacuum energy density, the effective cosmological constant that gives rise to 
dark energy. 

In fiat space-time, in the absence of galactic matter, the equation of motion is the non- 
linear Klein-Gordon equation (NLKG) 



(v^-^)0-A(l0r-Fo^)0 = O 



(6) 



with a conserved current density 

f = — (0*9^0 - 09^0*) = F^d^'a (7) 

In the presence of a galaxy (a generic term that includes a star), we add to the Lagrangian 
density an interaction term Ci^t, which represents the non-gravitational interaction between 
galaxy and scalar field. The galaxy is introduced as an external source, with a given energy 
current density J^. Its interaction with the scalar field is described phenomeno logically 
through a current-current interaction, as dictated by Lorentz covariance: 

.Cint = -VJ'J, (8) 

where r/ is a coupling constant, and J'^ is a four- vector: 

J" = (P, J) (9) 

where the given function p (x) describes the density profile of the galaxy. 
For a galaxy rotating as a rigid body with angular velocity 17, we have 

J = pQxr (10) 

where r is the distance vector from the center of the galaxy. In the equation of motion, 
1] and p occur only in the combination rjp. We define the coupling constant rj precisely by 
normalizing p : 

J,f.,^, (11) 

This coupling ultimately arises from particle-physics interactions, but we shall treat it phe- 
nomenologically by fitting rj to galactic data. (See later.) 
The NLKG now reads 

(v^ - ^) - A (101' - F^) - tvJ>^d,<p = (12) 

or 

(v'-^^<P-\{\<Pf-F^)<P-tVp(^^ + nxv.V<i^=0 (13) 

This equation contains a correlation length ,^o arising from the vacuum field Fq. 

«° - 7^ <"' 
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which we identify with the correlation length in galaxy clustering 15 1 

^0 ^ 11 Mpc = 3.3 X lO^^cm (15) 

This gives the mass of scalar-field quanta as 

mo = ^0^ ^ 10"^° eV (16) 

Another estimate of ^o can be obtained from the power spectrum of the CMB (cosmic 
microwave background). When analyzed in terms of Legendre polynomials Pi {cos 6), the 



angular distribution of thermal fiuctuations in the CMB shows a peak at £ ~ 200 [16||. Now, 
the angular pattern P^ (cos 6) divides the circle into 2i sectors, each subtending an angle 
A6 ^ n/i , and stuff within this angle should be correlated. So the correlation length would 
be of the order of ^o ~ 7ri?/^,where R ~ 10^^ cm is the distance of last scattering. Putting 
i = 200, we get ^o ~ 10^^ cm, which is not inconsistent with (HM . considering the qualitative 
nature of the argument. 

The NLKG here is a field equation, and is quite different from the relativistic N-body 
equations used in astrophysics. The scalar field here is the vacuum expectation value of a 
quantum field, which contains a high- momentum cutoff A. When A changes, the coupling 
A "renormalizes" , in order to maintain the identity of the theory. For the (p^ theory, A (A) oc 
In A. We take A = Aqcd + Acosmo, where Aqcd refers to the nuclear energy scale, and is 
taken to be time- independent, while Acosmo is the energy scale of the space-time metric, 
and decreases with the expanding universe. The latter is ignored because Acosmo << Aqcd- 
This is why we can treat A as a constant. 

III. GRAVITATIONAL INTERACTIONS 

We now include gravitational interactions in the Newtonian limit. The metric tensor is 
given by 



goo = -{l + 2U) (f/«l) 

^ 1 + 2U ^ ^ 

9jk = Sjk 

^ = VI + 2U (18) 
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where U is the gravitational potential. The equation of motion becomes 

-(l-2f/)|^ + V^0+^^ + Vf/-V0-A(|0r-Fo^)0-.r^p(^ + f^xr.V0)=O 

(19) 
The gravitational potential has contributions from the self-gravitation of the scalar field, 
and from the galaxy's gravitational attraction: 



\x — X 

nOO 



+ \^<l>\' + U\^\"-PoY 



dt 

Pgalaxy = ^galaxyP (20) 

where Mgaiaxy is the galaxy's total mass. With h = c = 1, the gravitational constant G is 
given by 

G= — (Planck length)^ = 3 x 10"^^ cm^ (21) 

An 

IV. GALACTIC DARK-MATTER HALO 

A galaxy immersed in the cosmic superfluid will become surrounded by a halo of higher 
superfluid density than the vacuum. This can be observed through gravitational lensing, 
and interpreted by us as dark matter. 

To see how the halo is formed, consider the situation in which a hypothetical galaxy had 
been formed in the vacuum scalar field long enough for the entire system to develop into a 
stationary state. We would have (p oc e^*"^*, with eigenfrequency u given by 

ur^riip) (22) 

where (p) denotes a spatial average of the galaxy's density. This formula can be derived in 



a variational formulation of the problem [17|]. The NLKG would then become 

V^^ - A ( 101^ - Fq^ - ^ j - TIP (cucp + inxr- V0) = (23) 

According to this equation, the vacuum field would be changed from Fq to Fi, with 

F? = F^ + ^ (24) 

9 



The situation described above is of course hypothetical. What actually happens is that a 
dark-matter halo will initially develop around the galaxy and spread, as illustrated schemat- 
ically in Fig.l. The initial speed of the spreading is the velocity of light, the character speed 
in the NLKG, but it will be slowed or stopped by gravitational interactions. To estimate 
the overall average speed, we assume that the halos around galaxies emerged with their 
formation, some 10^° years ago. We take the present radius to be 10 times galactic radius. 



or 10^ ly, judging from gravitational lensing pictures [l3|. This gives an average expansion 
speed less than 10~^c. 

The parameters Fi, u give rise to new correlation lengths: 

^' VXF, 

^. = - (25) 

CO 

and they are related to the vacuum correlation length through 

1 1 1 

72=72 + 71 (26) 

SI so Sw 

The region inside the dark- matter halo is governed by ^i, the outside region is governed by 
^0, and there is a transition region between the two. We may assume that the stationary 
equation fl23|) approximately holds within the halo, but of course it becomes invalid in the 
transition region and beyond. When the frequency u is sufficiently large, the halo region 
will become non-relativistic. 

For a disk-like galaxy, the halo would emerge as a surface layer on the disk, and grow in 
thickness under a force normal to the galactic disk, so as to lower the energy. Initially the 
layer would grow at the speed of light, the characteristic speed of the NLKG, but eventually 
it would be slowed or stopped by self gravitation of the layer, and the gravitational attraction 
from the galaxy. This could produce the "squashed beach ball" shape of the halo around 
the Milky Way y^]. The process is schematically illustrated in Fig.2. 

V. RELATION BETWEEN DARK ENERGY AND DARK MATTER 

The universe roughly consists of 3/4 dark energy and 1/4 dark matter. Since dark energy 
and dark matter are manifestations of the same scalar field, that fact imposes relations that 
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FIG. 1. Scalar-field profile showing the dark-matter halo around a galaxy. The vacuum scalar field 
has modulus Fq outside, and Fi > Fq inside. This creates a higher energy density in the halo, 
which can be observed via gravitational lensing. 
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FIG. 2. Formation of the dark-matter halo around a disk galaxy (side view). The halo initially 
grows at the speed of light normal to the galactic disk, and is eventually stabilized by gravity. 



can be tested. In the following we regard dark matter and halo as the same. Let 

Volume of all dark-matter halos 



Volume of universe 

Pdark matter 



(27) 



Pgalaxy 

where p denotes mass density. The following relation should hold between dark energy and 
dark matter, with neglect of observable matter: 



7Pdark matter + (1 — 7) Pdark 



energy 



Po 



(28) 
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where po is the energy density of the universe. The two terms above represent roughly 3/4 
and 1/4 of the total energy density of the universe, respectively [l9|. Thus 



^i "j j Pdark energy 
O'Pdark matter 
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37 . . 

Pdark energy Z Pdark matter '^ 'JTPdark matter V "/ 

1-7 

which leads to the relation 

J- Pdark energy ^ Pdark energy /on\ 

7 = o — = T7 (30) 

<J Pdark matter <J/ Pgalaxy 

Using the experimental values 

Pdark energy = 5.4 X 10*^ Cm""^ 

Pgalaxy = 1-8 X 10^^ Cm"^ (31) 

from ^] and [20| respectively, we get 

10"^ 
.^- (32) 

Since each galaxy occupies a volume the size of its halo, we should have 

7 = ^galaxy ( „ ^''^° ) (33) 

\ -''^universe / 

where -Runiverse is the radius of the universe, and A^gaiaxy is the number of galaxies. Equating 
the above to /^^10~^ and then solving for -Rhaio, we get 

J-, -fLunivcrsc -LU iU ly foA\ 

^halo = 1/3 .1/3 ~ TITT (34 j 

^^galaxy ■' ■' 

where we have used the experimental values -Runiverse = 5 x 10^^ cm, A^gaiaxy = 2 x 10^^ The 
average radius of a galaxy happens to be lO^ly, and the above can be rewritten -Rhaio ~ 

/ -fLgalaxy; OT 



f ^ ( ^^) (35) 

-Khalo / 



The definition ( 127|) is equivalent to / = (Mhaio /^galaxy) (-Rgaiaxy/-Rhaio) , where M denotes 



total mass. The above thus imp 
lensing, e.g. the "bullet cluster" 



ies that Mhaio/Mgaiaxy ~ 1- From images of gravitational 



13| . we estimate -Rhaio/-Rgaiaxy ~ 10, and obtain 

/ ~ 10-=^ (36) 



Thus, the halo is a dilute medium with mass density a thousand times smaller than that of 
the galaxy. 
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VI. SUPERFLUID VELOCITY AND QUANTIZED VORTICES 

A superfluid can flow pass walls without friction, as long as its velocity is below a critical 
velocity, according to the Landau criterion of superfluidity [5|. Similarly, a macroscopic 
object can move through a superfluid without friction, as long as its relative velocity is below 
the critical velocity, which is usually the phonon velocity associated with long-wavelength 
excitations, or Goldstone modes. For our cosmic superfluid, this is the velocity of light, and 
thus stars in a galaxy should be able to move frictionlessly through the dark-matter halo, 
since the latter is part of the superfluid. On the other hand, a superfluid can be induced to 
rotation about an axis, by creating quantized vortices in the superfluid (See below). 

In non-relativist ic theories, the superfluid velocity is deflned by Vg = (h/M) Va, where a 
is the phase of the complex order parameter, and M a mass. However, this is not bounded 
by the velocity of light, and clearly need modiflcation in the relativistic domain. In flat 
space-time, we generalize it as follows. We start with the conserved current j^ given by ([7]), 
which satisfles the continuity equation d^j'^ = 0, or 



dn 
dt 


+ V 


■J 


= 






j= 


=F^Va 






n 


= -F'& 



(37) 

where the fleld is represent in the form = Fe^'^, and a dot denotes time derivative. We 
interpret n as the superfluid density, and deflne the superfluid velocity Vg through 

f/Tl 

— + V ■ (nvs) = (38) 

This leads to the deflnition 



es = -- (39) 

a 

where ^g is a correlation length that is generally space-time dependent. Now the magnitude 
of Vs is bounded by the velocity of light, as long as d'^a is a time-like four vector. We recover 
the nonrelativistic case when a (r, t) = —Mt + v (r, t) , with M — )■ cxd. 
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It should be noted that Vs has physical meaning as a velocity field only when ^s is constant, 
or at least slowing varying. 

For a manifestly Lorentz-covariant formulation, we introduce a time-like unit four vector 

W^ = nd^a (40) 

with W^W^ = —1. Thus n is given by 

n = ^ (41) 

and 

w 



V2 

s 



W = -^=L= (42) 

To discuss vorticity, it is more convenient to work with d'^cr without the non-constant 
factor K, because it is d'^a that satisfies a quantization condition. Denote the spatial gradient 
by 

u = Va = Vs/^s (43) 

The line integral of u over a closed loop C must be a multiple of 27r, by continuity of the 
scalar field: 



ds-u= <b ds-Va = 2'kN (44) 

c Jc 



When A^ 7^ 0, there is a vortex line encircled by C. We refer the reader to Ref.[3| for a 
discussion of vortex dynamics. 

For disk-like galaxies rotating about an axis normal to the disk, vortex lines will be 
generated normal to the disk, and there is cylindrical symmetry at least near the surface of 
the galactic disk. The core size of the vortex line is proportional to the correlation length, 
which is ^1 inside the galactic halo, and increases as we go away from the center of the halo, 
and approach the larger vacuum value ^o outside the halo. The transition layer between 
the halo and the vacuum may be regarded as a fuzzy boundary where a vortex line can 
"terminate". This is illustrated in Fig.3. The halo of a galaxy with one vortex would look 
like a fat donut with a tiny hole. The picture might be different for a black hole, for which 
the strong space-time curvature may have important effects on vorticity. 
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FIG. 3. Core size of vortices created by a rotating galaxy are governed by the local correlation 
length. Inside the galactic it is some 30 times smaller than that in the vacuum. Consequently, the 
vortex cores increase in size towards the edge of the halo. 



We can roughly estimate the critical angular velocity for the creation of one vortex, using 
a cylindrical geometry. Consider a circular loop C of radius R, with an infinite vortex line 
normal to the loop through its center. The velocity u is tangential to the circle with constant 
magnitude uq. We find by explicit evaluation 



ds ■ U = 2'KRun 



(45) 



c 



Equating this to 2ttN yields uq = N/R, or 



R 



(46) 



The corresponding angular velocity is f] = Vs/R. Thus, to have one vortex inside radius R, 
the angular velocity must exceed the critical value 






(47) 



In Table 1, we give estimates of Qc for various stellar objects, putting ^s ~ R, i-e-j assuming 
that the correlation length of the superfiuid is greater than the size of the object. The actual 
angular velocities Q are also cited for reference. As we can see, a black hole could generate 
a million vortices, the Milky Way and neutron star are marginally critical cases, while the 
sun and earth fall short of criticality. 
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R (cm) ilc (rad/s) Q (rad/s) 

Milky Way lO^^ IQ-^^ iQ-is 

Earth 10^ 10 lO""^ 

Sun lO^i 10-1 10-"^ 

Neutron star 10^ 10"^ 10^ 

Black hole 10^ 10 10^ 

TABLE I. Critical angular velocity for creation of one vortex. Columns correspond to radius, 
critical angular velocity, actual angular velocity. 

When a black hole "drags" the surrounding superfluid into rotation, there would appear 
the order of 10^ vortex lines, which would form a cage around the black hole. This may 
explain the so-called "non-thermal filaments" observed in the neighborhood of Sagitarus 



near the center of the Milky Way |21 1 . 



VII. DETERMINATION OF MODEL PARAMETERS 

Assume that the correlation length ^i in a typical galactic halo is roughly the radius of 
the halo, which we take to be 10 times that of a typical galaxy: 

^1 ~ 5 X 10^2 ^^ (4g) 

This makes 

# = f^ ~ 30 (49) 

Since Fi >> Fq, we have Ff = F^ + u^/X^ ^ u'^/X, or 

CO ^ 1/6 (50) 

Now we calculate u using ( 122|) in the halo, taking F to be constant, and get 

UJ ^V{p) ^ Tj — (51) 

where Vhaio is the volume of the halo. Thus r] ~ wVhaio, or 

V ~ -^ ^6 (52) 
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The corresponding dimensionless parameter is 



^ = t4 = ^ 7^ -10 



a 3 V^o. 



2 



In out model, dark energy comes from Vq, and dark matter comes from deviations from 
the vacuum field: 



Pdark energy ^0 



Pdark matter 2 *^ 1 



37 (53) 



0, 

Since Fi >> Fq, we have 

XFf _ 1 

With dalD and dSOD, we get 

4 _ 2Fo _ o « ,, inl4„_-4 



(54) 



XF^ = ^l^ = 3.6x lO^^cm-4 (55) 
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We can find Fi and A separately by noting 

XI 
YFt 

The results are 



F^ = T^= (Ai^i^) e, (56) 



Fi ~ 10^° cm"^ 
Fo ~ 3.3 X 10^^ cm"^ 
A ~ 3.6 X 10-1°'^ (57) 

VIII. GOLDSTONE MODES AND VACUUM STABILITY 

Being a gravitating medium, the cosmic superfluid could be unstable against gravitational 
clustering. This would be manifested as the spontaneous emergence of patches of dark 



matter, analogous to the Jeans instability of a gravitating fluid J22|]. To investigate the 
stability of the vacuum, we perturbed it by exciting the Goldstone modes. The NLKG in 
the absence of matter, but with gravity in the Newtonian limit, is given by 

V^</>-(l-2f/)^ + ff + Vt/.V0-A(|0r-Fo^)0 = O (58) 
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where U is the gravitational potential satisfying the Poisson equation 

2 



V'U = AtiG 



+ |V0|V^ 



^o)' 



(59) 



The ground state corresponds to 0o = -^o- To investigate the Goldstone modes, we put 



Fo + f 



(60) 



where / is a small perturbation. When gravity is ignored, one can solve the linearized NLKG 
via the Bogoliubov transformation 



/ = ae-'^'"'-^-''^ + /3*e*('^*"''-^) 



(61) 



In this case, one obtains the dispersion relation w = fc, as dictated by Lorentz covariance. 
With inclusion of gravity, linearization in / leads to 



VV-(l-2t/)^ + — ^ + Vt/- V/-eo-' (/* + /) = 



(62) 



This is a nonlinear equation, because U depends on /. To solve it, we resort to an approxi- 
mation through the replacements 



dUdf _. df 
'dt'dt ^ ~'^'dt 

VU-Vf ^ -iq-Vf 



(63) 



where / on the right side is given by (EI]), and 7, q are real constants. We assume U « 1 
and neglect it, and obtain a set of homogeneous algebraic equations: 



Co-'a + [k' -00^ + ^0' + (q ■ k - ju)] (3 = 



(64) 



We average over the directions of q, and take (q ■ k) =0, ((q ■ k) ) 
determinant to zero, we obtain the dispersion relations 



\q^k'^. Setting the 



u^ = e + -(A± V2Bk^ + aA 



A = 2^,' + ^' 
B = q^ + 272 



(65) 
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FIG. 4. The Goldstone modes with uj = k are spht into "optical" and "acoustical" branches by 
gravity. In the acoustical branch, the sound velocity is pure imaginary in the region < /c < /cq, 
and this signifies vacuum instability, with formation of patches of dark matter of size Lq = 27r/A;o. 
In the present universe Lq is estimated to be much larger than the radius of the universe, and 
fco — ^ 0. Thus, the present vacuum should be stable against gravitational clustering. 



The Goldstone modes split into two branches (±) that are analogous to the "optical" (+) 
and acoustical (— ) modes in solids. As depicted in Fig.xx, both revert to the relation u = k 
when A; — )■ oo. At the low k end of the spectrum, u approaches a finite constant for for the 
optical mode, while it goes to zero at fc = fco for the acoustical mode, with the behavior 



u 



a/2A;o (k - ko 
{B - 2A) 



(66) 

li ko > 0, this signifies instability, in which the vacuum breaks up into patches of dark 
matter with characteristic size 



Lo 



27r 

kn 



(67) 



To estimate the instability length Lq, we use f l66|) in the form 



Lo = 2n{^-A 



-1/2 



-f-l 



(68) 



This is an implicit relation, because q depends on Lq. We can estimate q as follows. Suppose, 
in spontaneous fluctuation, the vacuum field develops a value Fi in a region of size Lq. The 
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energy density in this patch would be A {Ff — Fq), and the mass of the patch would be 

M = A {F,' - F,') Ll (69) 

By dimensional analysis, we estimate that 

q ~ f (70) 

Substituting this into fl68|) yields a quartic equation for Lq. The solution in the limit of weak 
gravity gives 

T O 

(71) 



^0 GF2 (Ft/F,' - 1) 
Using G = 3 X 10"^*^ cm^, Fq ~ 3.3 x lO^^ cm-\ we obtain 

^ = ^°" (72) 

We had obtained Fi/Fq ^ 30 in the galactic halo, but that was due to interaction between 
galaxy and vacuum field. Here, Fi is a property of the vacuum field, and is determined by 
balancing the weak gravity against the strong field potential. We expect that Fi/Fq ^ 1, 
which means Lq > 10^^ .^o- Since this length is overwhelmingly larger than the radius of the 
universe, we conclude that the present vacuum is stable against gravitational clustering. This 
of course does not preclude vacuum clustering during the early universe, when conditions 
were different. 

IX. NUMERICAL EXPLORATION 

We solve the NLKG with external-source galaxies whose densities have Gaussian distri- 
butions. The galaxies can rotate and move in prescribed ways, and our focus is on the 
response of the superfluid. Computations on moving galaxies are carried out in a 512 x 512 
two-dimensional spatial grid, with periodic boundary conditions. Those on vortex lattice 
creation are done in a 1024 x 1024 grid. For illustration purposes, model parameters are 
chosen arbitrarily. 

Fig. 5 shows a "film strip" of the superfluid in the presence of a moving galaxy, consisting 
of sequential plots of the scalar-field modulus. Frame 1 is the initial state, with velocity 
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directed upwards. Frame 2 shows the development of the dark-matter halo (white area) 
around the galaxy, and an outgoing transient cylindrical wave front in the superfluid. In 
frame 3, the transient wave continues to expand about the original center, while the halo 
moves with the galaxy like a soliton. In frame 4, the transient wave hits the boundaries. 

Fig. 6 shows the collision between two galaxies, initially made to move towards each other 
head-on. They continue on such a course, go through each other, and finally recede from 
each other. The frames show the evolution of the field modulus in contour plots, with the 
white areas corresponding to the dark-matter halos. The pictures show the merging of the 
two halos, the contortions of the combined halo, and finally the recession of the galaxies 
from each other, with their proper halos restored. 

Fig. 7 shows a double galaxy system rotating in a prescribed manner. In frame 2, two 
black dots appear, indicating the formation of vortices, and they continue to develop in the 
subsequent frames. It is the presence of these vortices that allow the superfluid to acquire 
angular velocity. 

Fig. 8 shows one moment during the development of a vortex lattice around a rotating 
galaxy. Such vortex lattices have been observed in cold-trapped atoms [;23|, and simulated 
in the NLSE J24|]. To our knowledge, this is the first simulation in the NLKG. The left 
panel is a contour plot of the field modulus, and the right frame is one of the phase of 
the field. The vortices form four rings, with the outer most ring located just outside the 
dark-matter halo (in white). The vortex counts of the rings are: 19, 12, 11, 9, a total of 
51, with the inner most ring being somewhat irregular because it is still being formed. In 
general, vortices are nucleated inside the galaxy, and migrate outwards. The phase contours 
show a radial pattern, in which the dark radial spokes represent "strings" across which the 
phase jumps by 27r, and these strings must be terminated by vortices. For example, the tips 
of the "chrysanthemum petals" coincide with the locations of the vortices on the outermost 
layer. Note that vortices in the outermost layer has much larger core size than those within 
the halo, because the correlation length is much larger. [See Fig. 3 and Eq. (H9l) .] 

Fig. 9 shows the number of vortices as function the angular velocity of the galaxy, for a 
range of the interaction parameter between galaxy and scalar field. 
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FIG. 5. Contour plots of the scalar-field modulus in 2D space, showing time evolution of the 
superfluid in response to the presence of a moving galaxy. The white areas correspond to the dark- 
matter halo. A transient cylindrical wave front can be seen propagating from the initial position 
of the galaxy, hitting the boundaries in the last frame. 




FIG. 6. Collision between two galaxies approaching each other along the vertical direction. The 
dark-matter halos merged, underwent deformations, and were finally restored when the galaxies 
recede from each other. 

X. THE NORMAL FLUID AND EXTENSIONS OF THE MODEL 



The NLKG presented here describes the vacuum scalar field interacting with an external 
source at absolute zero, and the field corresponds to a pure superfluid. At finite tempera- 
tures, a normal fluid component will emerge. At low temperatures, the normal fluid consists 
of a dilute gas of excitations that can be created from and absorbed back into the superfluid, 
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FIG. 7. A rotating double-galaxy system. The surrounding superfluid develops vortices (the black 
dots) in order to rotate with the galaxies. 

including 



phonons of the Goldstone modes described earlier; 

quanta of the scalar field, with mass given by (fT6|) . 

other particles that can be created from the scalar field: baryons, leptons, radiation, 
and possibly WIMPs (weakly- interacting massive particles) ll|, |l2| 



These particles would be described by classical distribution functions satisfying Boltzmann 
equations, which would be coupled to the NLKG in an appropriate manner. In such a 
description, the external source current J^ will be replaced by dynamical variables. 

Such investigations might throw light on extensions to the standard model. For illustra- 
tion, suppose the WIMP is a massive neutrino with mass generated by the vacuum scalar 
field, through a coupling of the form g (^iprcpip) , where r is some internal symmetry matrix. 
The mass of the WIMP would be qFq, where the vacuum field has the enormous value 
Fo ~ lO^^eV, according to ( 1571) . Thus, by dialing g, we can get a very wide range of masses. 

Even without considering the normal fluid component, our idealized system of NLKG 
with a perfect classical fluid as source can be generalized by endowing the classical perfect 
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FIG. 8. A vortex lattice surrounding a rotating galaxy. The left panel is a contour plot of the field 
modulus, showing four rings of vortices, the outermost of which lies beyond the galactic halo, and 
has a much larger core size. The right panel is a contour plot of the phase of the field. The radial 
spokes are "strings" across which the phase jumps by 27r. See text for a fuller description. 
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FIG. 9. Number of vortices around a galaxy, as a function of its angular velocity, for different 
values of the coupling constant between galaxy and scalar field. 



fluid with dynamics. One can add a dynamical part to the current: J^ = J^^ + q'^, where 
q'^ satisfles some Euler-like equation. This will address some immediate needs, for example, 
computation of galactic velocity curves. 



The earliest indication of dark matter came from the velocity curve of Adromeda 25| . 
which shows that "dust" surrounding the galaxy has a velocity much greater than that due 
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to the visible galactic mass. In the extension of the present model, one would be able to 
calculate the velocity, taking into account gravitational interactions with galaxy and halo, 
and coupling to the superfluid velocity field arising from possible vortices. Investigations 
along these lines are being pursued. 



Appendix A: Dimensionless form of NLKG 

For numerical computations, it is convenient to put the equation of motion into di- 
mensionless form. To do this, we scale all lengths with the vacuum correlation length 
Co = ( V^Fq ] , and introduce dimensionless quantities indicated with an overhead tilde: 

t = Cot 

X = C,oX 

n = ^^'n 

P = Co^'P 

V = ilv (Al) 



The equation of motion for the scalar field, with an external source p representing the density 
of a galaxy (or star), which is rotating with angular velocity f], is given in dimensionless 
form by 



<9V ^2 7 dUd^ 



90 



-(1 - 2U) ^+V <^+^^+V^-V0- ( "^ - M '^-'^P ( ^ + f^ X ? ■ V0 I = (A2) 

with the normalization convention 

d^xp = l (A3) 



The gravitational potential U (which is dimensionless) can be computed from Poisson's 
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equation: 



V^f/ = 4:7lG (pgalaxy + Ps. 
Pgalaxy -'"galaxyP 



Gps 



G 



dt 



V0 



+ 



(A4) 



The only adjustble paramete is the mass of the galaxy Mgaiaxy It is convenient to use the 
Milky Way as standard: 

GMMiikv wav = 1.2 X 10"^^ cm (A5) 



Let 



i Milky Way 



/3 



M 



galaxy 



Mt 



Milky Way 



then the Poisson equation takes the form 



^^U = iScip + C2 



dt 



+ 



V0 



2 1 

+ 2 



ci = 4.56 X 10"^° 
C2 = 10-1° 



(A6) 



(A7) 
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